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ABSTRACT 

HD 5980 is a unique system containing one massive star (star A) that is apparently entering 
the luminous blue variable phase, and an eclipsing companion (star B) that may have already 
evolved beyond this phase to become a Wolf-Rayet star. In this paper we present the results from 
FUSE observations obtained in 1999, 2000, and 2002 and one far-UV observation obtained by 
ORFEUS/'B^YS in 1993 shortly before the first eruption of HD 5980. The eight phase-resolved 
spectra obtained by FUSE in 2002 are analyzed in the context of a wind-eclipse model. This 
analysis shows that the wind of the emptor obeyed a very fast velocity law in 2002, which is 
consistent with the line-driving mechanism. Large amplitude line-profile variations on the orbital 
period are shown to be due to the eclipse of star B by the wind of star A, although the eclipse 
due to gas flowing in the direction of star B is absent. This can only be explained if the wind 
of star A is not spherically symmetric, or if the eclipsed line radiation is "filled-in" by emission 
originating from somewhere else in the system, e.g., in the wind- wind collision region. Except 
for a slightly lower wind speed, the ORFEUS/BEFS spectrum is very similar to the spectrum 
obtained by FUSE at the same orbital phase: there is no indication of the impending eruption. 
However, the trend for decreasing wind velocity suggests the occurrence of the "bi-stability" 
mechanism, which in turn implies that the restructuring of the circumbinary environment caused 
by the transition from "fast, rarefied wind" to "slow, dense wind" was observed as the eruptive 
event. The underlying mechanism responsible for the long-term decrease in wind velocity that 
precipitated this change remains an open issue. 

Subject headings: stars: individual (HD 5980) — stars: Wolf-Rayet — stars: binaries: eclipsing — stars: 
winds, outflows 
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1. Introduction 

HD5980 (Sk78, AV229) is the most Imninous 
eclipsing binary system in the Small Magellanic 
Cloud (SMC) and a member of the young stellar 
cluster NGC 346. It is an ideal object for studying 
interacting winds and the evolution of very mas- 
sive stars in binaries. The system achieved notori- 
ety in 1994, when its primary star (hereafter star 
A; M~50 M0) was discovered to be undergoing 
an eruptive event of a yet-undetermined nature. 
Although the system had been monitored since 
the early 1980s, the sudden brightening by ~2 
magnitudes in late 1993 went unrecorded except 
for the visual estimates of A. Jones (Bateson & 
Jones 1993). An intensive ultraviolet (UV) mon- 
itoring campaign was initiated in mid- 1994, and 
this database provides a unique record of the de- 
clining phase of the outburst. However, the most 
intriguing question that HD5980 raises involves 
the mechanism responsible for the eruption. We 
know that between 1987 and 1995 its UV spec- 
trum went successively through all the "late" ni- 
trogen Wolf-Rayet (WR) classes - WN6 WN7 

WN8 ^ WNll ^ WN6 - while also exhibit- 
ing significant changes in wind velocity and visual 
brightness; see Koenigsberger (2004) for a general 
review. Although some of the characteristics of 
the eruption are similar to those observed in Lu- 
minous Blue Variables (LBVs), both its extraor- 
dinary luminosity (^ 3 x 10^ Lq) and its increase 
in luminosity during eruption (~ 10'' Lq; Drissen 
et al. 2001) place it in a category currently shared 
only with rj Carinae. 

The second component of the eclipsing system 
(hereafter star B; M~28 Mq) also appears to be 
a WR star. It has been difficult to isolate its 
spectral characteristics from the combined spec- 
trum; see Koenigsberger (2004) for a discussion 
of this issue, which is further complicated by the 
possibility that some emission at line frequencies 
arises in a wind-wind colliding region (Moffat et al. 
1998). If star B is indeed a WNE star, then it is 
the evolved remnant of the originally more mas- 

^Postal Address: Dept. of Physics & Astronomy, The 
Johns Hopkins University, 3400 N. Charles Street, Balti- 
more, MD 21218. 

^ Based on observations made with the NASA-CNES- 
CSA Fax Ultraviolet Spectroscopic Explorer. FUSE is op- 
erated for NASA by the Johns Hopkins University under 
NASA Contract NAS5-32985. 



sive star of the binary system, and is therefore 
a good pre-supernova candidate. Alternately, it 
could simply be a less massive star whose outer 
layers have been significantly modified by a mass- 
transfer process, which has led to the current state 
of star A. 

The spectrum of HD5980 indicates the pres- 
ence of a third star. It can be seen in the con- 
tinuum (Breysacher & Perrier 1991), in the opti- 
cal spectra obtained in the 1980s (Niemela 1988), 
and in the "stationary" UV photosphcric lines that 
are visible in HST/STIS spectra (Koenigsberger 
et al. 2002). This component (hereafter star C) 
is likely to be an 04-6 supergiant that could be 
gravitationally bound to the star A-|-star B pair 
with a very long orbital period. Thus, it might 
be responsible for the instability of the system 
(Koenigsberger et al. 1994). Alternately, star C 
could simply be a line-of-sight object whose light 
contaminates the spectrum of the HD 5980 system. 
The presence of this third star complicates the 
interpretation of an already complicated system. 
However, because the evolution of massive stars 
plays such an important role in many astrophysi- 
cal phenomena, it is important to understand the 
mechanisms responsible for the peculiar behavior 
of HD 5980. 

In this paper, we analyze a set of phase-resolved 
spectroscopic observations of HD 5980 obtained 
with the Far Ultraviolet Spectroscopic Explorer 
{FUSE) satellite. These data were obtained over 
an interval spanning many orbital cycles, but in- 
clude a subset of 8 observations taken over ~6 or- 
bits, which permit orbital changes in the structure 
of the wind of star A to be assessed. In addition, 
we use an archival spectrum obtained by ORFEUS 
shortly before HD 5980 erupted. This fortuitous 
observation constrains the time scale for the devel- 
opment of the instability leading to the eruptions. 
Our analysis is organized as follows. In §2, we 
describe the observational material. The FUSE 
light curve is presented in §3, while the far-UV 
spectrum of HD 5980 and its orbital variability are 
discussed in §4. In §5 we analyze the orbital vari- 
ations in terms of a wind-eclipse model, and dis- 
cuss the implications. Additional discussion and 
our conclusions are presented in §6. 
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2. Observations and Data Analysis 
2.1. Observational Material 

HD 5980 has been observed several times with 
FUSE, though only those observations obtained 
during program P223 (PI: D. Massa) provide 
phase-resolved coverage. Program P223 produced 
8 observations that are well distributed in orbital 
phase, which were obtained over an interval of 
^116 days (i.e., ^6 consecutive orbital cycles) 
in 2002. Spectra from program P223 were sup- 
plemented by two snapshots of HD 5980, one of 
which was obtained as an early-release observa- 
tion (observation X0240202; PI: J. B. Hutchings), 
while the other was obtained as part of a PI- Team 
program to investigate the interstellar medium 
(observation P1030101; PI: K. R. Sembach). 

Table 1 provides details of the FUSE obser- 
vations of HD 5980. Successive columns list the 
name of the data set; the heliocentric Julian date 
(HJD) at the beginning of the integration; the to- 
tal integration time; the orbital phase at the start 
of the observation, which was computed from the 
ephemeris of Sterken & Breysacher (1997); the 
mean flux and the standard deviation in a fidu- 
cial continuum band; and the implied signal-to- 
noise ratio. All the FUSE spectra of HD 5980 
were obtained through the large (LWRS) aperture. 
Observation X0240202 was obtained in histogram 
(HIST) mode; all the others were obtained in time- 
tag (TTAG) mode. 

In addition to the FUSE spectra of HD5980, 
a spectrum obtained on 1993 September 17 by 
the Berkeley Extreme and Far-UV Spectrometer 
(BEFS) on-board ORFEUS was retrieved from the 
Multimission Archive at the Space Telescope Sci- 
ence Institute (MAST^). Details of this observa- 
tion arc also listed in Table 1. This spectrum 
was obtained only ~51 days before the eruption 
of HD 5980 began, and is the only UV spectrum 
that records the characteristics of the system at 
this crucial time. An additional stroke of luck is 
that the ORFEUS/BF,FS spectrum was obtained 
at (/) = 0.075, i.e., just after primary eclipse, when 
star A is "in front" of its companion. Hence 

^The archiving of non-HST data at MAST is supported by 
the NASA Office of Space Science via grant NAGS- 7584 
and by other grants and contracts. STScI is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 



the P Cygni absorption components in this spec- 
trum reflect the wind conditions of the emptor 
(star A). However, the ORFEUS/BEFS spectrum 
has poorer wavelength resolution and signal-to- 
noise (S/N) than the FUSE spectra. The spec- 
trum obtained from MAST was not processed fur- 
ther, although a wavelength shift of approximately 
-|-60 kms~^ was applied in order to align interstel- 
lar features with their counterparts in the FUSE 
spectra. 

Finally, as described below, our analysis of 
HD5980 was supplemented by FUSE spectra 
of two other stars in the SMC: Sk80 (AV232; 
07 laf-h; see Crowther et al. 2002) and Skl08 
(AV332; WN4:-h06.5 I:; see Mallouris 2003). 
These observations were obtained in TTAG mode 
through the LWRS (Sk80) or MDRS (Skl08) 
aperture. Additional details for these observa- 
tions are also provided in Table 1. 

2.2. Processing of FUSE Spectra 

The instruments on-board the FUSE satellite 
consist of four aligned, prime-focus telescopes and 
Rowland circle spectrographs (Moos et al. 2000; 
Sahnow et al. 2000). Two of these "channels" 
arc coated with LiF to provide high efficiency 
longward of ~1000 A, while the other two are 
fabricated from SiC to give high throughput be- 
tween the Lyman limit and ^1100 A. The spec- 
tra are recorded by two photon-counting detectors, 
each of which consists of two segments. Thus, a 
FUSE observation produces 8 independent spec- 
tra, which are labelled by their channel (LiFl, 
LiF2, SiCl, SiC2) and their detector segment (A 
or B). Except for a small region near 1100 A, the 
entire FUSE waveband is redundantly covered by 
2-4 of these spectra. 

The FUSE spectra were uniformly extracted 
and calibrated with the standard reduction pipeline 
(CalFUSE version 2.4.1). Subsequent processing 
steps combined the 8 individual spectra from a 
given observation into a single spectrum. Since the 
different spectra have inherently different quality 
(due to the different effective areas of each chan- 
nel) and suffer from instrument artifacts to differ- 
ing degrees (e.g., thermally induced channel mis- 
alignment), these manipulations were conducted 
interactively in three steps. 

First, residual zero-point shifts in the wave- 
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length scales for the different spectra were re- 
moved by comparing the observed positions of 
H2 lines formed in the interstellar medium of the 
Galaxy with the laboratory wavelengths listed by 
Jenkins & Peimbert (1997). The corrections de- 
rived in this way were applied to each spectrum 
individually. The corrections for data from the 
SiC channels were always smaller than 0.25 A, 
while for LiF data they were in general smaller 
than 0.03 A. The overlapping regions of spec- 
tra from different channels and detector segments 
were checked carefully to ensure that the resultant 
wavelength scales were consistent. When neces- 
sary, the spectra were trimmed to minimize the 
effect of nonlinearities in the wavelength scale due 
to residual distortion near the edges of detector 
segments. 

Second, the aligned spectra were inspected vi- 
sually. In general, only the highest quality data for 
any particular wavelength interval were retained 
for the merged spectrum. Gaps in wavelength cov- 
erage were filled with data of poorer quality only 
when necessary. LiF2A spectra were favored over 
LiFlB spectra for coverage longward of ~1100 A 
because the fiux calibration of LiF IB spectra is 
severely compromised by an optical artifact known 
as "the worm" (Sahnow 2003). 

Third, the aligned, trimmed, and hand-selected 
spectra from the various channels and detector 
segments were combined by using the IRAF'^ rou- 
tine scombine. The spectra contributing to the 
final, combined data product were assigned equal 
weights (w — 1) except for SiClA spectra, which 
were given w=0.25 because of their poorer qual- 
ity. Finally, the data were rebinned to a constant 
wavelength step of 0.05 A. 

Seven of the eight P223 observations of HD 5980 

were combined to produce a global mean refer- 
ence spectrum for the P223 data set. This aver- 
age spectrum has S/N ~40 per data sample in the 
continuum near 1000 A. Observation P2230103 
(0 =0.124) was excluded because detector 1 was 
not operating at full voltage for the entire expo- 
sure, so that spectra extracted from it (i.e., LiFlA, 
LiFlB, SiClA, SiClB) were anomalously noisy. 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 



2.3. Far-UV Continuum Fluxes 

Continuum placement in spectra of WR stars 
is very uncertain for all UV wavelengths due to 
the number and breadth of emission and P Gygni 
wind profiles. Shortward of ~1000 A, the problem 
is aggravated by the large number of interstellar 
absorption lines, particularly from H2, and by the 
confluence of the Lyman series of atomic H. In ad- 
dition, the spectrum of HD 5980 contains absorp- 
tion from a velocity system associated with the 
foreground supernova remnant SNR0057— 7226.^ 

Despite these complications, we identified sev- 
eral spectral bands that allow the continuum level 
to be defined rehably: e.g., 1010-1020 A and 
1055-1060 A. In particular, the wavelength inter- 
val 1059.6-1060.6 A is free from interstellar and 
photospheric lines. We used the average value 
of the flux over this interval, which we denote 
by ^Aioeo, to characterize continuum variations of 
HD 5980 in the far-UV. The values of Fxweo, along 
with the root-mean-square scatter over the contin- 
uum band [a-(FAio6o)], and the corresponding S/N 
ratio are provided in Table 1. 

3. Far-Ultraviolet Light Curve of HD 5980 

The orbital period of HD 5980 is 19.2654±0.0002 
(Sterken & Breysacher 1997). As indicated in Fig- 
ure 1 (taken from Koenigsberger 2004), orbital 
phase = 0.0 is defined by the minimum of the 
primary eclipse; i.e., when star A eclipses star B. 
Owing to the substantial eccentricity of the or- 
bit (e = 0.32,0.27,0.30 according to Breysacher 
& Perrier 1991; Moffat et al. 1998; Kaufer et al. 
2002, respectively), the secondary eclipse occurs 
at (/) = 0.36. The radius of star A prior to erup- 
tion (Ra(1978) = 21 Rq) and the radius of star B 
(Rb = 15 Rq) were determined by Breysacher & 
Perrier (1991) from analysis of the visual light 
curve. Since the orbital inclination is ^ 86° 
(Breysacher & Perrier 1991), the visual contin- 
uum eclipse at </> =0.0 is total. However, the 
current radius of star A is uncertain. In 1994 De- 
cember, Koenigsberger et al. (1998a) determined 

"'In the 1200-1800 A spectral range, the velocity system 
of the SNR is at Vheiio = +312 ± 3 kms^^ and +343 ± 
3 kms^^, and includes lines for ions ranging from O I to 
C IV (Koenigsberger et al. 2001). The same system has 
been detected in FUSE spectra of the O VI resonance dou- 
blet (Hoopes et al. 2001; Danforth et al. 2003). 
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it to be Ra(1994) = 35 Rq. An analysis of the 
contemporary light curve is needed to confirm the 
expectation that the radius of star A has been de- 
creasing systematically over the intervening years. 

Our FUSE observations permit a far-UV light 
curve to be derived, which covers both eclipses. 
The top panel of Figure 2 plots i^Aioeo from Ta- 
ble 1 as a function of orbital phase. The curve 
superposed on the data illustrates the predicted 
eclipse depths from the model described in §5. 
The bottom panel of Fig. 2 shows the differen- 
tial magnitude of HD 5980 with respect to Sk 108 
(large triangles), and also indicates the visual light 
curve determined by Breysacher & Perrier (1980, 
crosses). Although the number of far-UV data 
points is limited, there is generally good agreement 
between the far-UV and optical light curves, ex- 
cept for the FUSE data obtained in 1999 and 2000 
(indicated by parentheses). Evidently the system 
was systematically brighter in these early observa- 
tions compared with 2002, presumably because it 
was still declining from the maximum brightness 
attained during its 1994 eruption. The depths of 
the optical and far-UV eclipses at (f> = 0.36 are 
similar, while at (/> = 0.00 the far-UV eclipse seems 
to be deeper than the optical eclipse. However, 
given the large optical photometric variations de- 
tected on timescales of ~6 hours, little more can 
be concluded from this single point at (j) —0. 

The ORFEUS/BEFS measurement of the far- 
UV continuum flux is the point in Fig. 2 at 
(j> = 0.075. Although it has the largest uncer- 
tainty of any of the i^Aioeo measurements, its 
mean value agrees very well with that obtained 
at nearly the same orbital phase from FUSE spec- 
trum P2230102. This agreement of pre- and post- 
outburst measurements separated by 8.9 years 
(168 orbital cycles) is very interesting, since it 
suggests that HD 5980 has returned to its pre- 
out burst state. In this context, it is interest- 
ing to note that even though the ORFEUS/BEFS 
spectrum was obtained shortly before the start of 
the eruption, the FUV continuum fluxes did not 
presage dramatic changes. 

4. Far-UV Spectrum and Variability 

The average FUSE spectrum of HD 5980 is illus- 
trated in Figures 3 and 4, where it is also compared 
with the spectrum of Sk 108 (dotted hues). Sk 108 



is an ideal template for the analysis of HD 5980 
because: 1) it is also in the SMC and is thus 
similarly affected by reddening; and 2) it consists 
of a WN4 star and an 06.5-supergiant (FocUmi 
et al. 2003), and should resemble the combined 
spectrum of star B -|- star C. A striking difference 
between these objects is the number and inten- 
sity of emission lines in the far-UV spectrum of 
HD 5980. While the FUSE spectrum of Sk 108 has 
strong emission only in the O VI resonance dou- 
blet, the spectrum of HD 5980 exhibits emission 
features due to S VI, N IV, N III, He II, O VI, 
S IV, P V, and C IV. 

Assuming that the combined spectrum of star 
B -|- star C resembles the spectrum of Sk 108, then 
the source of the strong emission lines in HD 5980 
must be star A and, possibly, the wind-wind col- 
lision (WWC) region. At longer UV wavelengths 
(1200 - 1700 A); most emission lines arc also dom- 
inated by star A (Koenigsberger 2004), although 
star B appears to contribute significantly to the 
P Cygni features at N V A1240, C IV A1550, and 
He II A1640. One feature that is associated ex- 
clusively with star A (and, possibly, the WWC 
region) is the semi-forbidden N IV] A1486 line, 
which displays extremely large line-profile varia- 
tions. These variations are usually interpreted in 
terms of a non-spherically symmetric line-emitting 
region. 

In the FUSE wavelength region, the P V 
AA1117, 1128 doublet presents a similar oppor- 
tunity for isolating emission arising only in star 
A, with the advantage that these resonance lines 
produce strong P Cygni absorption components. 
For 0-type stars, the intensity of the P V doublet 
is correlated with the stellar luminosity (Walborn 
et al. 2002), and the great strength of these lines 
in the spectrum of HD 5980 is consistent with the 
large luminosity derived for star A (Koenigsberger 
et al. (1998b, La ~ 3 x IO^Lq)). 

We used the standard deviation of 7 of the 8 
spectra obtained in 2002 to quantify the spec- 
tral variability of HD 5980. As shown in Figs. 3 
and 4, the largest values of this dispersion coin- 
cide with the positions of strong emission lines. 
At most other wavelengths, the time variations 
have a base level of 0.5 x lO"""^^ erg/s/cm^/A, 
which refiects the continuum variations due to 
the eclipses but also includes cycle-to-cycle fluc- 
tuations. Smaller variations occur only at wave- 
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lengths that coincide with saturated absorption 
lines; e.g., the strong absorptions due to interstel- 
lar Lyman /3 at ~1025 A. 

The P V All 17 line is only moderately af- 
fected by blending with other stellar and inter- 
stellar lines, and can thus be used to study line- 
profile variability in more detail. Figure 5 illus- 
trates the orbital variations of P V A1117 for 
the eight FUSE spectra obtained in 2002. In 
this montage, the spectrum obtained near primary 
eclipse (P2230101; (j) = 0.001) is used as a tem- 
plate against which variations in the other profiles 
are compared. In order to make this comparison 
meaningful, two systematic effects were allowed 
for. First, the profiles were shifted in velocity so 
that they represent the emission in the frame of 
reference of star A. The radial velocity curve of 
Niemela et al. (1997) was used to remove the or- 
bital motion of star A. Second, the profiles were 
re-normalized with respect to the sum of the com- 
bined continua of star A -|- star B. This was ac- 
complished by removing the light of star C, which 
was assumed to contribute approximately 33% of 
the total continuum in the far-UV, as it does in 
the visual region of the spectrum (Breysacher & 
Perrier 1980). By taking the observed continuum 
level at phases 0.124 and 0.775 as representative 
of the sum of all three continua, the contribution 
from star C was estimated to be ~ 2.1 x 10~^^ 
erg/s/cm^/A in the vicinity of the P V doublet. 
This value was subtracted from each of the spectra 
in turn. Finally, the 8 P V profiles were renormal- 
ized by dividing by the out-of-eclipse continuum 
level of star A + star B (4.2 x 10~^^ erg/s/cm^/A). 

Fig. 5 indicates that the primary cause of line- 
profile variability in P V is the appearance of 
strongly enhanced absorption (with respect to 
(f> = 0.001) in the P Cygni profile when star A 
is "in front" of star B. Conversely, when star A is 
eclipsed by star B at w 0.360, the low-velocity 
part of the P Cygni absorption trough is weakened 
compared with primary eclipse. This clear pattern 
of behavior can be attributed to wind and stellar 
eclipses, as demonstrated by the model discussed 
in the next section. 



5. Predicted Line Profile Variations from 

Eclipses 

P Cygni profiles in the spectrum of a hot star 
are produced in an extended and rapidly expand- 
ing stellar wind. Variability in these lines can be 
used to diagnose: a) the presence of eclipses; b) 
aspherical geometry of the wind due, e.g., to col- 
lisions with the outfiow from a companion star; 
and c) structural changes in the density, velocity, 
or ionization conditions in the wind. Eclipses and 
wind-wind collisions produce periodic variability 
on orbital timescales. Wind eclipses, in partic- 
ular, may be used to derive information about 
the geometry and structure of the stellar winds 
of the components of an early-type binary system 
(Miinch 1950; Willis et al. 1979; Koenigsberger 
1990; Auer & Koenigsberger 1994). 

Figure 6 provides a representation of the star A 
-|- star B binary system at four orbital phases, un- 
der the assumption that both stars possess stellar 
winds. The geometry of the shock cone produced 
by the collision of these outflows is approximated 
with the formalism developed by Canto et al. 
(1996), with stellar wind velocities estimated by 
Koenigsberger (2004), t;oo(star A) = 2000 kms"\ 
Woo (star B) = 2600 kms~^ and estimates for the 
mass-loss rates of Ma = 1 x 10"'' Mq/jt, Mb = 
2 X 10~^ M0/yr. These mass-loss rates should 
be taken as representative values only, since the 
only reliable determination was made shortly after 
the eruption in 1994 (Koenigsberger et al. 1998b; 
Drissen et al. 2001) and may be larger than dur- 
ing quiescent times. Furthermore, the wind veloc- 
ity has varied significantly since the 1980s. We 
assume that the wind of star A consists of an 
inner region extending out to a distance Vacceh 
within which the atmosphere accelerates and ap- 
proaches the maximum wind speed, Voo- In Fig- 6, 
Taccei = 1.4Ra is represented by the circle traced 
by long dashes. Beyond VacceU the wind expands 
at the constant speed v^o- The wind that flows in 
the direction of the companion undergoes a highly 
supersonic shock when it encounters the compan- 
ion's outfiow, which brakes the radial component 
of its motion and forces it to move along the work- 
ing surface of the shocks. The kinematic proper- 
ties of the shock surface are determined by the 
conservation of momentum in the hydrodynami- 
cal system; see, e.g.. Canto et al. (1996). 
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Fig. 6 illustrates that continuum radiation from 
star A or star B can traverse several physically dis- 
tinct regions before leaving the system. Consider, 
e.g., the elongation illustrated in Fig. 6(b). Start- 
ing from star B, the short, dashed lines parallel to 
the line-of-sight cross in succession the wind of star 
B, the shock cone surface, and finally the wind of 
star A before reaching a distant observer located 
at the bottom of the figure. All these regions have 
velocity components in the direction of the ob- 
server, and therefore all contribute to the shape 
of the line profile at negative velocities, i.e., to the 
morphology of the absorption trough of a P Cygni 
profile from the binary system. Hence, while the 
emission component of the P Cygni profile con- 
tains contributions from large volumes of wind ma- 
terial, the contribution to the absorption trough is 
confined to the wind contained in the columns of 
material projected onto the stellar disks. 

The discussion of wind-profile variability can 
be simplified if we choose a P Cygni line that is 
formed in the wind of star A but not in the wind of 
star B. We argued above that the WN4:+06.5 I: 
binary system Skl08 should be a good template 
for the combined spectrum of star B and star C, 
since star B is believed to be a WNE star. Sk 108 
lacks P Cygni emission due to P V and thus we 
assume that the entire P Cygni feature of P V in 
HD 5980 is attributable to star A. In this case, we 
expect to observe the intrinsic P Cygni absorption 
component of star A at orbital phase (j) =0.00, 
when it fully eclipses its companion. At orbital 
phases immediately before and after (p =0.00, the 
continuum radiation of star B must traverse re- 
gions of the wind of star A to reach the observer. 
This is the "wind eclipse". Because the wind of 
star A has velocity components in the direction 
of star B as well as in the direction of the ob- 
server, the wind eclipse causes an enhancement of 
the P Cygni absorption component and a reduc- 
tion in the intensity of the emission component. 
The latter is due to absorption at positive veloci- 
ties from the wind that is expanding toward star 
B. At 4> =0.36, when star B eclipses star A, the 
P Cygni absorption component should be signifi- 
cantly reduced, since the column of wind material 
from star A where this absorption is produced is 
occulted by star B. 



5.1. Model Calculation 

These concepts can be quantified by c;omput- 
ing synthetic line profiles with a simple radiative 
transfer code that includes the effects of contin- 
uum and wind eclipses (Auer & Koenigsberger 
1994; Flores et al. 2001). This program treats 
the radiative transfer using the Sobolev approx- 
imation only where this approximation is valid, 
and performs the exact integration of the radia- 
tive transfer equation in wind regions where the 
velocity gradient is too small to satisfy the Sobolev 
criterion. For simplicity, the model assumes spher- 
ically symmetric line-forming regions and neglects 
both wind-wind collision effects and scattering 
processes. The model profiles thus provide base- 
line information on the eclipse effects alone. Dis- 
crepancies between the observed and predicted 
line profile variations can be used to assess the 
importance of other phenomena not incorporated 
in the model (e.g., the wind of star B; wind- wind 
collisions; or an aspherical wind geometry). 

As before, we assume that the wind obeys a 
linear velocity law out to a distance Vacceh beyond 
which the velocity has the constant value t^oo-^ In 
addition, we assume that the ionization structure 
is such that P V is present to a distance r^ax- 
The computer program allows for the possibility 
of a radial dependence of the ionization struc- 
ture through the parameter /(r), which is given 
by f{r) = x('')^^('')^^) where x(r) is the opacity. 
Through a trial-and-error process, the free param- 
eters of the model were adjusted until the general 
features of the P V A1117 line profile at (j) =0.001 
were reproduced. Emphasis was placed on trying 
to match the shape of the P Cygni absorption com- 
ponent. This was achieved with the set of input 
parameters labeled "MpOl-1" in Table 2. Figure 7 
compares the synthetic line profile computed with 
these parameters for (/> = 0.001 (solid line) with 
the observed profile. 

A satisfactory match between the synthetic and 
observed P V P Cygni profile of star A could only 
be achieved with a very fast wind velocity law. 
Table 2 shows that the terminal speed, Vmax, of 
model MpOl-1 is achieved within raced = 1.4Ra. 
The consequences of varying raced are illustrated 

^ All distances are expressed in units of the radius of star A, 
. All velocities are in units of the "turbulent velocity" , 

'^turb- 
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in the profiles plotted at the bottom of Fig. 7: a 
slower velocity (corresponding to raccei = 1-6 Ra) 
produces a P Cygni absorption trough that is too 
deep at low velocities (short dashes), while a faster 
velocity law (corresponding to raccei = 1-3 Ra) 
produces too little low-velocity absorption (long 
dashes). Since rapid acceleration (i.e., a "fast" 
velocity law) is a key feature of winds driven by 
radiation pressure in spectral lines, these results 
suggest that the current wind of star A is driven 
by this mechanism. 

For a fixed opacity structure, the maximum ex- 
tent of the line-emitting region, Tmax^ determines 
the intensity of the line profile, since a larger vol- 
ume contains a larger number of line-emitting ions. 
The consequences of varying rmax are illustrated 
by the profiles in the upper section of Fig. 7. 
Since the long- and short-dashed curves bracket 
the plausible range of fits to the observed profile, 
we conclude that 3.3 Ra < Tmax < 4.5 Ra for the 
region emitting P V A1117 radiation. 

The value of the terminal wind speed, Voo= 
1750 kms"-'^, constrains the product {vmax) x 
{vturb)- However, small values of Vturb produce 
very sharp and deep absorption at Voo, contrary 
to what is observed. Instead, we find that val- 
ues of Vturb of ~70-80 kms^^ (which corresponds 
to ~5% of Woo) are required to match the soft 
blue edge of the absorption trough. As a re- 
sult, the "edge" velocity (i.e., the velocity where 
the blue edge of the absorption trough meets the 
continuum) is significantly faster than ?;oo- This 
is a characteristic of strong P Cygni absorption 
troughs in 0-type stars in general (Prinja et al. 
1990). The occurrence of soft blue absorption 
edges is generally attributed to the presence of 
additional velocity dispersion in the wind, which 
is frequently characterized as "turbulence" . How- 
ever, despite its widespread occurrence, the origin 
of this velocity dispersion is not known. 

The other significant discrepancy between the 
model and the observed line profiles is the sharp 
emission peak near line center; see Fig. 7. The 
red component of the P V doublet at 1127 A pro- 
vides a built-in check of this discrepancy. Since 
the doublet separation (+2693 kms~^) is greater 
than Voo for star A (1660 kms""'^; see Table 3), 
the two components are not radiatively coupled 
and can be treated independently. The oscilla- 
tor strength of red component of the P V doublet 



is half that of the blue component, so the model 
profile for P V A1127 was calculated with half the 
opacity, /(r), listed for MpOl-1 in Table 2, but 
holding all other parameters fixed. Figure 8 shows 
the synthetic profile for the doublet obtained by 
adding the computations for each component to- 
gether. It confirms that the synthetic profiles pro- 
vide a good representation of the observed line 
strengths, and that the sharp peak near line center 
is the most discrepant feature for both components 
of the doublet. (Note that the strong absorption 
at ~1000 kms~^ in the absorption trough of the 
red component is an interstellar feature.) We will 
show below that this problem may be solved if we 
assume that star C is responsible for the excess 
emission near velocity. 

The model parameters listed in Table 2 also 
permit line profiles to be computed for different 
orbital phases. In the model, phases arc character- 
ized in terms of the impact parameter, p, which is 
defined as the minimum distance between the line- 
of-sight to star B from the center of star A. For 
illustrative purposes, p can be read directly from 
the abscissae in Fig. 1. The phases for the model 
computations were chosen to have impact param- 
eters corresponding to the phases of the FUSE oh- 
servations. 

The phase-resolved progression of synthetic line 
profiles is shown in Figure 9, where each profile is 
labeled with the corresponding impact parameter 
(expressed in units of Ra)- The line profile com- 
puted for primary eclipse {p — 0) is used as tem- 
plate, to show the predicted range over which ex- 
cess absorption should be present. The left panel 
corresponds to orbital phases when star B is oc- 
culted by star A, while the right panel shows the 
phases when the portions of the wind of star A 
are eclipsed by star B. The profiles in both panels 
can be compared directly with their counterparts 
in Fig. 5. 

Fig. 9 shows that the P Cygni absorption 
trough of P V in the spectrum of star A should 
decrease in strength near secondary eclipse, while 
the emission lobe is relatively unaffected. Thus, 
the line profile behavior shown in Fig. 5 is ex- 
plained qualitatively by the combined effects of 
wind and physical eclipses. There is, however, 
one important feature of the theoretical predic- 
tion that does not appear in the observations: for 
small impact parameters (e.g., p » 2Ra), the 
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model predicts enhanced absorption over the en- 
tire velocity range from —Vao to +^00, projected 
onto star B. In contrast, the observations only 
show excess absorption from — Wqc to w « 0. 

Thus, it appears that either: (a) the wind of 
star A flowing toward star B does not have the 
same structure as the wind emerging from the op- 
posite side of star A; or (b) the predicted excess 
absorption on the "red" side of the line profile 
is filled in by emission from some process that 
is not currently included in the model. For (b), 
this emission must arise from gas that is not pro- 
jected against either star A or star B, and which 
is receding from the observer at orbital phases 
-^O.OibO.lS. A plausible candidate that would 
produce emission is the WWC region illustrated 
in Fig. 6(a). However, the opening angle of the 
shock cone would need to be sufficiently small for 
the projection of the WWC flow velocities along 
line-of-sight to be as large as +1000 kms~^. In 
addition, this putative emission would have to be 
visible at the opposite orbital phases ~ 0.36) 
at similar speeds approaching the observer. Al- 
though the observed proflle at this orbital phase 
does display what could be interpreted as excess 
emission, the same feature can also be interpreted 
in terms of a physical eclipse of the absorption- 
forming column of material in the wind of star A. 
This effect is predicted by the model: see, e.g., 
the p = 0.3 Ra proflle in the right-hand panel of 
Fig. 9. Thus, it seems more likely that the wind 
emitted by opposing hemispheres of star A has 
substantially different structure. 

5.2. Correcting for Line Emission From 
star C 

If star C is a mid-0 type supergiant (Koenigs- 
berger ct al. 2002), its spectrum should contain 
prominent P Cygni profiles in the resonance lines 
of P V; see, e.g., Walborn et al. (2002). We used 
an archival FUSE spectrum of Sk80 (07 Iaf+) as 
a template for star C, and subtracted it from each 
spectrum of HD 5980 in order to compensate for 
the contamination it introduces. Since the con- 
tinuum level of SkSO matches the estimated con- 
tinuum strength of star C (i.e., ~33% of the total 
recorded for HD 5980), this operation also corrects 
the FUSE spectra for the contribution to the con- 
tinuum from star C. Figure 10 illustrates the 8 line 
profiles from FUSE program P223 after subtract- 



ing the spectrum of Sk80 and correcting for the 
orbital motion of star A (as in Fig. 5). 

A new set of models was computed for the 
P V A1117 profile of star A at = 0.001. Match- 
ing parameters to achieve a good fit with the ob- 
served profile yielded once again raccei = 1-4 Ra, 
but required larger values of r^ax and f{r). The 
adopted input parameters are listed in Table 2 un- 
der the label "MpOl-2". The increase in rmax and 
/(r) compensate for the fact that the P Cygni 
profile of the P V resonance line in the spectrum 
of Sk80 has a prominent, low- velocity absorption 
component. Subtracting this profile from the spec- 
trum of HD 5980 produces more emission, and con- 
sequently the model requires more optical depth 
over a larger volume to achieve a good match. 

Figure 11 shows the P V A1117 profile of star A 
at (/) = 0.001 after the removal of the contribution 
from star C. The "best- fit" profile computed with 
the parameters of model MpOl-2 is shown as a 
continuous line, which corresponds to rmax = 5-3. 
For comparison, models with rmax = 5.6 Ra and 
^max = 4.8 Ra are also shown. The match is sig- 
nificantly better than the one illustrated in Fig. 7, 
especially in the emission lobe near v = 0. Thus 
we conclude that this feature is indeed an artifact 
of stationary "third light" contamination by star 
C. 

It is important to note that the assumed P V 
contribution from star C is overestimated by 
adopting Sk 80 as a template. This line's emission 
is maximum in supergiants at 06—07 (Walborn 
et al. 2002), while star C is likely to be 04-06. 
Thus, star A's actual wind structure is between 
that given by models MpOl-1 and Mp02-2, al- 
though the uncertainty does not affect the con- 
clusions regarding the wind accelerating region. 
Sk 80 is the "earliest" 0-supcrgiant in the SMC 
for which there exists a FUV spectrum, and thus 
provides the only template we can use at this time. 

The phase-dependent line profile variability 
predicted by this second set of calculations is 
shown in Figure 12. The variations are quali- 
tatively similar to those generated from model 
MpOl-1; compare, e.g., with Fig. 9. The main dif- 
ference is the more pronounced wind eclipse due 
to the larger physical extent of the P V emitting 
region which, in this case, extends beyond the 
orbital radius. 
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5.3. Absolute Dimensions of the System 

A theoretical light curve can be constnicted 
from the model continuum intensities that were re- 
quired to achieve a good match with the observed 
line profile of P V A1117. An additional, unknown 
requirement is the size of star B relative to star A. 
By varying this single parameter until a good fit 
with the far-UV light curve (Fig. 2) was achieved, 
we determined that Kb = 0.7 Ra in 2002. By fur- 
ther adopting Rb =15 R© (Breysacher & Perrier 
1991), we infer that the far-UV continuum of star 
A was Ra =21 Rq in 2002. With this value of 
Ra, the length scale for all model parameters is 
established, which in turn allows the distances in 
Figures 1 and 6 to be calibrated in units of R©. 
Table 3 lists the derived parameters for the wind 
of star A in 2002. 

It is interesting to note that the value of R^ 
that we find from the FUSE light curve obtained in 
2002 is the same as that deduced by Breysacher & 
Perrier (1991) from optical data obtained in 1978. 
This coincidence suggests that the optically thick 
radius marking the base of the stellar wind of star 
A has not changed, despite the disruptions caused 
by the eruptive phenomenon that took place in the 
system. 

6. Discussion and Conclusions 
6.1. The Wind Structure of Star A 

We have used phase-resolved wind profiles of 
the P V resonance doublet obtained by FUSE in 
2002 together with a simple wind eclipse model to 
constrain the velocity law of star A, the relative 
continuum intensities of star A and star B, and 
the physical dimensions of the HD 5980 system. 
We find that, to a first approximation, the P V 
wind profiles require that the wind of star A is 
expanding with a "fast" velocity law, as is typical 
of outflows driven by radiation pressure in spectral 
lines. Systematic variations of profile strength and 
morphology as a function of orbital phase are well 
reproduced by a combination of "wind eclipses" of 
star B and physical eclipses of star A by star B. 

However, there are four important discrepan- 
cies between the observations and the synthetic 
line profiles obtained from our simplified model. 
The first consists of the presence of an emission 
"spike" near the rest velocity in the components 



of the P V resonance doublet that is not predicted 
by the model. We show in §5.2 that this problem 
can be resolved by assuming that the excess emis- 
sion arises from the third component of the sys- 
tem, star C. The presence of a P V wind feature 
is consistent with the classification of star C as an 
04 06 supergiant on the basis of its photospheric 
absorption lines (Koenigsberger et al. 2001). Un- 
fortunately, solutions for the three remaining dis- 
crepancies are less straightforward. 
Absence of a wind eclipse at t; > : The wind 
eclipse around phase ~0.0 produces reduced emis- 
sion over the range of velocities associated with 
the column of material projected against the disk 
of star B. We have shown that the wind of star 
A expands rapidly, so that it achieves its terminal 
speed before it reaches the radius of its relative 
orbit about star B. Thus, at phases immediately 
before and after (/) =0.00, the entire emission lobe 
should be significantly reduced in intensity by the 
wind eclipse. However, the observations show the 
expected effect only at w < kms~^; i.e., only on 
the blue side of the P Cygni profile. 

Thus, it seems as if only the portion of the wind 
of star A that is approaching the observer has the 
geometry and kinematics assumed by the model. 
We see two possible explanations for this discrep- 
ancy: 

1 . There is an extended region of P V emission 
associated with star B that compensates for 
the absorption produced by the wind eclipse. 
This emission could come from the surface 
of the shock cone associated with the WW C 
region, which effectively truncates the spher- 
ical symmetry of the outflow. 

2. The wind of star A is not spherically sym- 
metric, but is significantly perturbed in the 
direction of star B. However, the asymme- 
try must be present in the wind accelera- 
tion region, which we have shown lies very 
close to the surface of star A, far from the 
expected location of the WWC stagnation 
point. Hence, if this explanation is correct, 
we conclude that the wind structure of star 
A in the direction of its companion differs 
from that of other directions, as opposed to 
a wind structure that is truncated by the 
WWC region. This difference could either 
be intrinsic to star A or could be the result 
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of "sudden radiative braking" (Gayley et al. 
1997). Note that if the wind structure to- 
wards the companion is non-standard, then 
the WWC surfaces that arc drawn in Fig. 6 
are not vahd. For example, a slower wind to- 
wards the companion would have the effect 
of moving the stagnation point towards star 
A. 

Strength of P Cygni Absorption at (f) =0.36: 
If P V A1117 is formed exclusively in the wind of 
star A, then the eclipse model predicts that the 

absorption trough of its P Cygni profile should 
exhibit significant weakening during secondary 
eclipse, because star B occults the column of ma- 
terial responsible for producing this absorption. 
However, only a modest degree of weakening is 
observed at slow velocities; see, e.g.. Fig 5 and 
Fig. 10. At large velocities, v > v^o, the absorp- 
tion seems to be even more pronounced during 
secondary eclipse. A natural explanation for this 
behavior is that the wind of star B also emits and 
absorbs P V; i.e., that our assumption concerning 
the unique association of P V with star A is not 
correct. Unfortunately, the present version of the 
wind-eclipse model is not able to test this idea, 
since a rigorous calculation must include detailed 
radiative transfer through two stellar winds and 
the WWC region. Instead, we simply note that 
if the spectrum of star B does posses a P Cygni 
wind profile in the P V resonance doublet, then it 
differs significantly from the spectrum of the WR 
star in Sk 108. 

Extent of the P Cygni absorption edge: The 

final discrepancy is the extent of the P Cygni ab- 
sorption edge. Although Uoo= 1750 kms~^ for 
star A, the edge velocity is v^dge > —2200 kms^-'^. 
The greater extent of Vedge is a well-known phe- 
nomenon in 0-type and WR stars, so it is not too 
surprising to find it in HD 5980 as well. What is 
surprising, however, is that the extended, "soft" 
blue edge of the P Cygni absorption trough does 
not vary with orbital phase. Since star B also 
drives a strong stellar wind, it is possible that 
its radiation field further accelerates the outflow 
from star A that has already achieved Voo- How- 
ever, since this speculation cannot be tested di- 
rectly with the present data set, the issue must 
remain open for now. 



6.2. New Insights into the Eruption Mech- 
anism 

The changes in the wind of HD5980 prior to 
the eruptions in 1993-1994 were characterized 
by a progressive decrease in terminal velocity, 
which was accompanied by a systematic increase 
in density (Koenigsberger et al. 1998b). The first 
sign that a significant change was occurring in 
HD 5980 can be found in its lUE spectra of 1986 
(Koenigsberger 2004), although peculiarities were 
already present in the early 1980s (Niemela 1988). 
The amplitude of the perturbation seems to have 
grown gradually over the subsequent ~13 years 
until some kind of critical state was reached. This 
critical state produced the sudden eruption in 
1993, which was followed by a second, stronger 
eruption in 1994. According to the visual data of 
Albert Jones (reproduced in Koenigsberger 2004), 
the eruption started around HJD 2,449,299 (1993 
November 7), which is only 51 days (2.6 orbital 
cycles) after the ORFEUS/BEFS spectrum was 
acquired. 

Figure 13 compares the P V A1117 wind profile 
in the ORFEUS/BEFS spectrum with its counter- 
part in the FUSE P2230102 spectrum. Both spec- 
tra were obtained at essentially the same orbital 
phase, though they are separated by 168 orbital 
periods. In addition to slightly stronger P Cygni 
absorption around —500 kms^^, the only mor- 
phological difference between the two spectra is 
the location of the discrete absorption component 
that marks the position of terminal speed: in 1993 
Voo= 1530 kms-\ while Voo= 1750 kms'^ in 2002. 
Hence, the trend for decreasing wind velocity per- 
sisted from the early 1980s up to -^51 days before 
the start of the eruption. 

Evidently, the overall structure of the wind of 
star A was very similar just prior to the first erup- 
tion in 1993 and in 2002, when HD 5980 was in 
a quiescent state. Hence, we infer that the wind 
must have been driven by the same mechanism 
at both epochs; and we have argued on the ba- 
sis of the rapid acceleration required to model the 
P V wind profile that in 2002, this mechanism 
was radiation pressure in spectral lines. Conse- 
quently, in 1993 the conditions in the star must 
have changed from supporting a "normal" radia- 
tivcly driven outflow to a state of significantly en- 
hanced mass loss in less that 51 days; i.e., what- 
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ever set of critical conditions initiated the eruption 
occurred within a remarkably short time interval. 
This rapid reconfiguration of the properties of the 
wind provides a strong constraint on the unknown 
mechanism responsible for the eruption. 

The fact that the wind velocity decreased sys- 
tematically for many years before the eruption 
suggests that some critical limit was reached in 
the autumn of 1993. If the decreasing wind veloc- 
ity is associated with an increasing wind density, 
as appears to be the case from the growth of emis- 
sion line intensities over the 1980-2000 timescale 
(Koenigsberger 2004), it is possible that the criti- 
cal conditions are related to the "bi-stability limit" 
(Lamers et al. 1995; Vink et al. 1999). In this 
case, the eruptions would be interpreted as the ob- 
servable manifestation of the transition from the 
"fast, rarefied wind" side of the bi-stability limit 
to the "slow, dense wind" side and the concomi- 
tant re-structuring of the circumbinary environ- 
ment. Since two major eruptions occurred within 
~1.5 years, it appears that star A remained very 
close to its "bi-stability limit" for this long before 
beginning to relax to its pre-outburst state. 

Of course, the crucial unanswered question is 
what underlying evolutionary process caused star 
A to move towards its "bi-stability limit" in the 
first place, e.g., by triggering enhanced mass- loss 
rate or decreasing the vc^Iocity of the radiatively 
driven outflow. Identifying this trigger mechanism 
in HD 5980 is certain to improve our understand- 
ing of eruptive phenomena that occur in other 
massive stars and binary systems. Since these out- 
bursts may drive significantly greater mass loss 
than possible via stellar winds, they have great 
potential to alter both the evolutionary histories 
of these stars and the yields of chemically enriched 
material they provide to their local environments. 

GK thanks Albert Jones for providing his visual 

magnitude estimates. This investigation was sup- 
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and by CONACYT grant 36569. 

REFERENCES 

Auer, L. H., & Koenigsberger, G. 1994, ApJ, 436, 
859 

Bateson, F. M., & Jones, A. F. 1993, Pub. R. As- 
tron. Soc. New Zealand, 19, 50 



Breysacher, J., & Perrier, C. 1980, A&A, 90, 207 

— 1991, in lAU Symp. 143, Wolf-Rayet Stars 
and Interrelations with Other Massive Stars in 
Galaxies, edited by K. A. van der Hucht, & 

B. Hidayat (Dordrecht: Kluwer), 229 

Canto, J., Raga, A. C., & Wilkin, F. P. 1996, ApJ, 
469, 729 

Crowther, P. A., Hillier, D. J., Evans, C. J., Fuller- 
ton, A. W., De Marco, O., & Willis, A. J. 2002, 
ApJ, 579, 774 

Danforth, C. W., Sankrit, R., Blair, W. P., Howk, 
J. C., & Chu, Y.-H. 2003, ApJ, 586, 1179 

Drissen, L., Crowther, P. A., Smith, L. J., Robert, 

C, Roy, J.-R., & Hillier, D. J. 2001, ApJ, 546, 
484 

Flores, A., Auer, L. H., Koenigsberger, G., & Car- 
dona, O. 2001, ApJ, 563, 341 

Foellmi, C, Moffat, A. F. J., & Guerrero, M. A. 
2003, MNRAS, 338, 360 

Gayley, K. G., Owocki, S. P., & Cranmer, S. R. 

1997, ApJ, 475, 786 

Hoopes, C. G., Sembach, K. R., Howk, J. C, & 
Blair, W. P. 2001, ApJ, 558, L35 

Jenkins, E. B., & Peimbert, A. 1997, ApJ, 477, 
265 

Kaufer, A., Schmid, H. M., Schweickhardt, J., & 
Tubbesing, S. 2002, in ASP Conf. Ser. 260, In- 
teracting Winds from Massive Stars, edited by 
A. F. J. Moffat, & N. St-Louis (San Francisco: 

ASP), 489 

Koenigsberger, G. 1990, A&A, 235, 282 

— 2004, Revista Mexicana de Astronomia y As- 

trofisica, 40, 107 

Koenigsberger, G., Auer, L. H., Georgiev, L., & 
Guinan, E. 1998a, ApJ, 496, 934 

Koenigsberger, G., Georgiev, L., Peimbert, M., 
Walborn, N. R., Barba, R., Niemela, V. S., 
Morrell, N., Tsvetanov, Z., & Schulte-Ladbeck, 
R. 2001, AJ, 121, 267 

Koenigsberger, G., Kurucz, R. L., & Georgiev, L. 
2002, ApJ, 581, 598 



12 



Koenigsberger, G., Moffat, A. F. J., St-Louis, N., 
Auer, L. H., Drissen, L., & Seggewiss, W. 1994, 
ApJ, 436, 301 

Koenigsberger, G., Pena, M., Schmutz, W., & Ay- 
ala, S. 1998b, ApJ, 499, 889 

Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, 
D. M. 1995, ApJ, 455, 269 

Mallouris, C. 2003, ApJS, 147, 265 

Moffat, A. F. J., Marchenko, S. V., Bartzakos, 
P., Niemela, V. S., Cerruti, M. A., Magalhaes, 
A. M., Balona, L., St-Louis, N., Seggewiss, W., 
& Lamontagne, R. 1998, ApJ, 497, 896 

Moos, H. W., et al. 2000, ApJ, 538, LI 

Miinch, G. 1950, ApJ, 112, 266 

Niemela, V. S. 1988, in ASP Conf. Ser. 1:, Progress 
and Opportunities in Southern Hemisphere Op- 
tical Astronomy, edited by V. M. Blanco, & 
M. M. Phillips (San Francisco: ASP), 381 

Niemela, V. S., Barba, R. H., Morrell, N. L, & 
Corti, M. 1997, in ASP Conf. Ser. 120, Lumi- 
nous Blue Variables: Massive Stars in Transi- 
tion, edited by A. Nota, & H. Lamers, (San 
Francisco: ASP), 222 

Prinja, R. K., Barlow, M. J., & Howarth, L D. 
1990, ApJ, 361, 607 

Sahnow, D. J. 2003, Proc. SPIE, 4854, 610 

Sahnow, D. J., et al. 2000, ApJ, 538, L7 

Sterken, C., & Breysacher, J. 1997, A&A, 328, 269 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 
1999, A&A, 350, 181 

Walborn, N. R., FuUcrton, A. W., Crowther, P. A., 
Bianchi, L., Hutchings, J. B., Pellerin, A., Son- 
neborn, G., & Willis, A. J. 2002, ApJS, 141, 
443 

Willis, A. J., Wilson, R., Macchetto, F., Bccck- 
mans, F., van der Hucht, K. A., & Stickland, 
D. J. 1979, in The First Year of lUE, edited by 
A. J. Willis (London: University College Lon- 
don), 394 



This 2-column preprint was prepared with the AAS lAT^TjX 
macros v5.2. 



1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 t 1 1 1 1 1 1 1 1 1 t 1 1 1 1 1 

e = 0,27. cj=133'- 




0,47 J 








\ 0.25 


■ 0.66,^ 


13 0.20 


:o,74|, 


W 0.16 


\ 

\^ 

: 0,82\ 


G) - 

V_V 1 0.13 

/o.io J 


'. 0.89° 


\ /^.07 '. 




0,97 0,01 

. . . , 1 . . . . 1 . . . . 1 . . . . 1 , . . . 1 . . ■ 



-160 -100 -50 50 100 150 

distance (Rg) 

Fig, 1, — Schematic representation (from Koenigs- 
berger 2004) of the orbit of star A in a reference 
frame that is fixed with respect to star B, viewed 
from above the orbital plane. The observer views 
the system from a perspective that is nearly in the 
plane of the page, looking "up" from the bottom 
along a line through impact parameter (abscissa) 
0. The open squares are located at 20° intervals 
and are labelled by their orbital phases. 
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Fig. 2. — Far-UV continuum light curve of 
HD5980. Top: the average continuum flux _F\i06o, 
expressed in units of 10^^^ erg/s/cm^/A and 
not corrected for reddening. Error bars indicate 
c(Fxio6o); see Table 1. The solid curve corre- 
sponds to the prediction of the model described 
in §5. Bottom; the magnitude difference between 
HD 5980 and Sk 108 (triangles) compared with the 
visual light curve in 1978 (crosses; Breysacher 
& Perrier 1980). F measurements from 1999 
and 2000 are enclosed in parentheses. Evidentally 
HD 5980 was systematically brighter in the far-UV 
at these epochs. 
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Fig. 3. — The short-wavelength region of the mean 
far-UV spectrum of HD5980 in 2002. The light 
tracing indicates the mean spectrum of HD 5980; 
see §2.2. The dark lines show the standard devi- 
ation about the mean spectrum caused by orbital 
variations. This "variation spectrum" has been ar- 
bitrarily shifted by -2 x 10"-^^ erg/s/cm^/A. The 
dotted hue shows the FUSE spectrum of the SMC 
binary Skl08 (WN4:-h06.5 L), while a normal- 
ized absorption-line spectrum from a Kurucz LTE, 
line-blanketed model atmosphere is plotted at the 
bottom (shifted by-2.5 x 10"^^ erg/s/cm^/A). 
Major features are identified. 
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Fig. 4. — Same as Fig. 3 for the long-wavelength 
region of the FUSE waveband. The emission fea- 
ture marked by an "X" is due to scattered solar 
light. 
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Fig. 6. — Geometry of the orbit and WWC region 
for different phases of the HD 5980 system. The 
observer is at the bottom in the four panels. Ver- 
tical dashes indicate the columns of wind material 
that produce P Cygni absorption troughs, which 
includes regions where enhanced absorption is pro- 
duced due to wind eclipses. The circular broken 
lines around star A show the extent of the accel- 
erating wind (1.4R^). The dotted circles show 
r = 4.1 Ryi, which marks the estimated lower limit 
of the P V line emitting region; see §5. 
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Fig. 5.— Line profiles of P V A1117 at differ- 
ent orbital phases, compared to the line profile at 
(j) =0.001 (dotted Hues). The wind eclipse (left 
panel) only affects approaching velocities. 
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Fig. 7.— The observed line profile of P V A1117 
at (j> = 0.001 (dotted) compared with the profile 
synthesized from model MpOl-1 (light solid). Dis- 
continuous lines show computed line profiles with 
all but one model parameter held constant. Top: 
Effect of varying the maximum extent of the line- 
emitting region. Short dashes: r^ax = 3.3 Ka', 
long dashes: r^ax ~ 4.5 Ka- Bottom: Effect 
of varying the wind velocity law. Short dashes: 
Taccei = 1.6 Ra; long dashcs: raced = 1.3 Ra- 
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Fig. 8. — Comparison between observed and syn- 
thesized profiles of the P V AA1117, 1128 doublet 
at (j) =0.001. Excess emission is apparent near the 
center of both lines. 
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Fig. 9. — Orbital line-profile variations computed 
from model MpOl-1. Individual profiles are la- 
beled the value of the impact parameter (ex- 
pressed in units of R^) of the line-of-sight to star 
B from the center of star A, and correspond to the 
orbital phases covered by the FUSE observations 
(Fig. 5). 
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Fig. 10.— Line profiles of P V A1117 at differ- 
ent orbital phases, compared profile at =0.001 
(dotted lines). In contrast to Fig. 5, the spectrum 
of Sk80 has been subtracted to account for the 
contributions of star C. The sharp emission excess 
seen near w = in Fig. 5 is no longer present. 
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Fig. 11.— Model fit to the P V All 17 profile after 
subtraction of the spectrum of Sk80, which was 
used as a proxy for the contributions of star C. 
The continuous line is the closest match to the 
data, which was computed with the parameters 
of model MpOl-2 (Table 2). The long and short 
dashes correspond to the same model but with 
Tmax = 5.6 Ra and Vmax = 4.8 Ra, respectively. 
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Fig. 13. — Comparison between P Cygni profiles 
of the P V All 17 resonance line observed by OR- 
FEUS/BEFS in 1993 (dashed hue; (j) = 0.075) and 
FUSE in 2002 (soUd line; (f> = 0.084). The OR- 
FEUS/BEFS spectrum was shifted in velocity by 
—65 kms^^ to align nearby interstellar lines and 
by +0.5 X 10^^^ erg/s/cm^/A to align the contin- 
uum levels. The velocity of the discrete absorption 
component at the time of the ORFEUS observa- 
tion was 220 kms~^ slower than at the time of the 
FUSE observation, but otherwise the profiles are 
very similar. 



Fig. 12. — Orbital line-profile variations computed 
from model MpOl-2. 
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Table 1 
Journal of Observations 



Observation HJD*^ Int. Time [ks] Phase'' Faio6o'' cr(i?Aio6o)'^ S/N 



HD 5980 




















BEFS1068'= 


49248, 


.0358 


1, 


.810 


0, 


.075 


5.73: 


0.74: 


7 


X0240202 


51472, 


.3964 


3, 


.153 


0, 


.535 


7.33 


0.43 


17 


P1030101 


51728, 


.0736 


5, 


.745 


0, 


.806 


6.61^ 


0.31 


21 


P2230105 


52432, 


.2891 


4, 


.342 


0, 


.360 


4.84 


0.37 


13 


P2230107 


52478, 


.8300 


3, 


.953 


0, 


.775 


6.16'' 


0.39 


16 


P2230108 


52481, 


.7160 


8, 


.050 


0, 


.925 


6.02 


0.28 


22 


P2230101 


52483, 


.1789 


6, 


.957 


0, 


.001 


4.45 


0.26 


17 


P2230102 


52484, 


.7828 


6, 


.925 


0, 


.084 


5.75 


0.27 


21 


P2230103 


52485, 


.5381 


7, 


,894 


0, 


.124 








P2230104 


52488, 


.8000 


7, 


.871 


0, 


.293 


5.97 


0.30 


20 


P2230106 


52548, 


.5025 


3, 


.886 


0, 


.392 


5.03 


0.36 


14 


Sk80 




















P1030201 


51728, 


.2801 


11, 


.700 






2.25 


0.10 


23 


Skl08 




















P1030401 


51728, 


.5619 


13, 


.958 






3.82 


0.14 


27 



''^HJD - 2,400,000.0 at the start of the integration. 

'^Orbital phase computed from the ephcmcris of Stcrkcn & Breysacher (1997). 

"^Mean continuum flux between 1059.6 and 1060.6 A, except for Sk 80 where the 
interval 1057—1058 A was used, in units of 10~^^ erg/s/cm^/A. 

Standard deviation in the mean continuum flux. 

•^Observation with ORFEUS/BEFS. 

^There is a broad dip at 1058.54 A at this phase. 
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Table 2 

Parameters for P V All 17 Model Calculation 



Parameter 


MpOl-1 


MpOl-2 


Description 




4.1 


5.3 


Maximum extent of emitting region 


Voo 


25.0 


20.8 


Terminal speed of the wind 


'^turb 


70 


80 


"Turbulent" speed in kms~^ 


'^accel 


1.4 


1.4 


Radius over wliieli the wind is accelerating 


fo 


400 


400 


Opacity factor at the stellar surface 


f{r),ri 


2600, 1.3 


4100, 1.3 


Opacity factor out to ri 


f{r),r2 


220, 4.1 


220, 5.3 


Opacity factor out to r2 




48, 32 


48, 32 


Continuum intensity of each star 


'^orb 


4.76 


4.76 


Mean orbital separation 


Rb 


0.7 


0.7 


Radius of star B in units of R^ 



Table 3 

Wind Parameters of Star A in 2002 



Parameter 


Value 


Comments 


Rmax 


111 Rq 


Maximum extent of emitting region = rmax Vturb 




1660 kms"^ 


Terminal speed of the wind of star A, =VaoVturb 


'^turb 


80 kms-i 


"Turbulent" speed in the wind in star A 


^accel 


29 Rq 


Radius of accelerating wind region for star A, = TacceiRA 


Ra 


21 Rq 


Radius of star A in 2002 
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